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Abstract 

In a previous paper we showed that, for any n > m + 2, most sets of n points in are determined 
(up to rotations, reflections, translations and relabeling of the points) by the distribution of their pairwise 
distances. But there are some exceptional point configurations which are not reconstructible from the 
distribution of distances in the above sense. In this paper, we present a reconstructibility test with running 
time 0{n^^). The cases of orientation preserving rigid motions (rotations and translations) and scalings 
are also discussed. 

Introduction 

In this paper, we present a quick and easy (but slightly imperfect) solution to the problem of characterizing 
the shape of sets of n points in Euclidean space, so-called n-point configurations, for any positive integer n. 
More precisely, an n-point configuration is a collection of n points in M'". Point configurations often arise in 
biological and medical imagery, as well as in the fields of archaeology, astronomy and cartography, to name 
just a few. For example, stellar constellations, minutiae of fingerprints, and distinguished points (landmarks) 
on medical images represent point configurations. 

An important problem of computer vision is that of recognizing point configurations. In other words, 
the problem is to determine whether two point configurations have the same shape, that is to say, whether 
there exists a rotation and a translation (sometimes a reflection and/or a scaling are allowed as well) which 
maps the first point configuration onto the second. Let us first concentrate on the case of rigid motions, 
i.e. rotations, translations and reflections in R™. Note that any rigid motion can be written as (A/, T), where 
M is an orthogonal m-by-m matrix and T is an m-dimensional (column) vector 

One of the biggest difficulties in trying to identify point configurations up to rigid motions is the absence 
of labels for the points: one does not know, a priori, which point is going to be mapped to which. If the points 
were already labeled in correspondence, then, following the so-called Procrustes approach (Gower [7]), one 
could analytically determine a rigid motion which maps the first string as close as possible (in the sense, 
for example) to the second. The statistical analysis of such methods is presented in Goodall [6]. Another 
way to proceed would be to compare the pairwise (labeled) distances between the points of each point 
configurations (Blumenthal [2]). Indeed, the following well known fact holds. See, for example, Boutin and 
Kemper [4] for a simple proof. 

Proposition 0.1. Let pi, . . . ,p„ and qi, . . . ,qn be points in M™. If\\pi~pj\\ ~ \\<li~ Qj \\ for every — 
1, . . . , 71, then there exists a rigid motion (M, T) such that Mpi + T = qi, for every i = 1, . . . , n. 

A variety of methods have been developed for labeling the points of two n-point configuration in corre- 
spondence. See, for example. Hartley and Zisserman [8] for a description of some of these methods. But 
labeling the points is a complex task which we would much rather do without. Invariant theory suggests a 
possible approach for recognizing unlabeled points. The idea consists in comparing certain functions of the 
pairwise distances between the points of the configuration which have the property that they are unchanged 
by a relabeling of the points. These are often called graph invariants and have been computed in the case of 
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Value 


# of Occurrences 


1 


4 


V2 


2 



Table 1 : Distribution of distances of a unit square. 

n = 4 by Aslaksen et al. [1], and n = 5 by the second author [5, page 220]. Unfortunately, the case n — 6 
or larger still stands as a computational challenge. Moreover, the invariants used are polynomial functions of 
the distances whose number and degrees increase dramatically with n. They are thus very sensitive to round 
off errors and noise. 

In the following, we study an alternative approach based on the use of a very simple object: the distri- 
bution of the pairwise distances. The distribution of the pairwise distances of an n-point configuration is an 
array which lists all the different values of the pairwise distances between the points in increasing order and 
the number of times each value occurs. For example, the distribution of distances of four points situated at 
the corners of a unit square is given in Table 1 . 

Obviously, such a distribution remains unchanged under any rigid motion of the point configuration as 
well as any relabeling of the points. For n = 1, 2 or 3, it is easy to see that the distribution of distances 
completely characterize the n-point configuration up to a rigid motion. For n > to + 2, we proved that, 
most of the time, this distribution completely characterizes the shape of the point configuration (see [4, 
Theorem 2.6]). 

To simplify our discussion, we introduce the concept of reconstructibility from distances. 

Definition 0.2. We say that the n-point configuration represented by pi, . . . ,pn G M™ is reconstructible 
from distances if, for every . . . , G having the same distribution of distances, there exists a rigid 
motion (M, T) and a permutation n of the labels {1, . . . , n} such that Mpi + T = qT^^i), for every i = 
l,...,n. 

In the following, we shall often identify a point configuration and one of its representation pi , . . . , p„ e 
M"*. This is done for simplicity and we hope it will not create any confusion. Please note that the question 
of reconstructing the point configuration from its distribution of distances will not be addressed in this paper 
We suspect this is quite a challenging problem. In fact, our guess is that this problem lies in the complexity 
class NP; it might even be NP-complete. 

Theorem 2.6 of Boutin and Kemper [4] actually implies that there exists an open and dense subset ft C 
(M™)" of reconstructible point configurations. In Section[2 we concentrate on the planar case m = 2 and 
give an algorithm in 0{n^^) steps to determine whether a point lies in O. (A simple Matlab implementation 
of this algorithm is given in the appendix.) A generalization to other dimensions to is also mentioned. 
Section 3 describes how an additional distribution can be used in the planar case in order to compare the 
orientation of two point configurations. In Section 4, we show that a slightly modified distribution can be 
used to completely characterize most point configurations up to rigid motions and scalings. 

1 Reconstructible configurations 

Denote by V the set of pairs 

'P = ^ J,i, j = 1, • • • ,n} 

Consider the group of permutations 5^^^ of the elements of V. For any E S^n^^ and any G V, we 

denote by (p ■ {i, j} the image of {i, j} under ip. For two point configurations to have the same distribution 
of distances means that there exists a permutation ip E Sf^n-^ which maps the labeled pairwise distances of 
the first configuration onto the labeled pairwise distances of the second configuration. More precisely, if 
pi, . . . ,Pn e and qi, . . . ,qn £ R™ have the same distribution of distances, let dij — \\pi — and 
d'i j = \\qi ^ QjW^^ for ^11 £ ^- Then there exists ip e Sj-n-j such that 

'^vlhj} = 't*'-^J' ^ 



RecoDstructible Point Configurations 



3 



For close enough point configurations, we have proved in [4] that one does not need to keep track of the 
labehng of the points. The proof is very short and we reproduce it here for completeness. 

Proposition 1.1. For any n-point configuration pi, . . . ,pn € M™, there exists a neighborhood U of 
(pi, . . . ,Pn) (z (M™ ) " such that i/ (gi , . . . , g„ ) G U is an n-point configuration with the same distribution 
of distances as that of {pi , . . . , p„), then the two point configurations are the same up to a rigid motion and 
a relabeling of the points. 

Proof. Let us assume the contrary. Then there exists a sequence of n-point configurations {q\, . . . , qn}kLi 
converging to pi, ... ,p„, and a sequence of permutations {(ySfcj^j^ C S'^n-j such that none of the q\,...^q^ 
can be mapped to pi, . . . ,p„ by a rotation and a translation and a relabeling, but the distances — 
\\pi - pjW^ are mapped to the distances ^j, = \\q^ - by ipk so d^^.{i^j} = d^'- .^j, for all {i, j} e V. 
By taking a subsequence, we may assume that tpf^ = is the same for every k since 5^™^ is a finite group. 

Taking the limit, we obtain that d^.^ijy — limfc_^oo for E V. By continuity of the distance, 

this implies that d^.^i_jj = d^^^y, for all £ V. Therefore, ~ d'^- -y for every G V. By 

Proposition lO.il this implies that q^, . . . ,q^^ and pi, . . . ,pn are the same up to a rigid motion, for every k, 
which contradicts our hypothesis, and the conclusion follows. □ 

Unfortunately, the size of the neighborhood is unknown and varies with the points pi, . . . ,Pn, so this 
local result is not very practical. We now consider the global case. Observe that some of the permutations in 
S^n-j correspond to a relabeling of the points. More precisely, ip corresponds to a relabeling of the points if 

there exists a permutation tt : {1, .... n} ^ {1, . . . , n} of the indices such that i^- {i, j} — {7r(i), 7r(j)}, for 
every {i,j} € V. Relabelings are the good permutations: if the permutation mapping the labeled pairwise 
distances of a point configuration onto the labeled pairwise distances of another configuration is a relabeling, 
then the two configurations are the same up to a rigid motion. We need to know what distinguishes the good 
permutations from the bad permutations. The following lemma, which is central to our argument, says that, 
informally speaking, a permutation is a relabeling if it preserves adjacency. 

Lemma 1.2. Suppose n ^ A. A permutation (p g SfJ^•^ is a relabeling if and only if for all pairwise distinct 
indices i,j,k£ {1, . . . , n} we have 



Proof. For n < 3, every Lp e 5^™^ is a relabeling, and the condition M.W is always satisfied. Thus we may 
assume n > 5. It is also clear that every relabeling satisfies ( II. H . 

Suppose that ip e S^n-^ is a permutation of V which satisfies ( II. H . Take any k, I G {1, . . . , n} 

pairwise distinct and assume, by way of contradiction, that p ■ {i, j} D (p ■ {i, k} O (p ■ — 9. Then the 
injectivity of ip and the condition ( ILH imply that we can write ip ■ {i,j} = {a, b}, p ■ {i,k} = {a, c}, and 
(p ■ {i, 1} = {&, c} with a,b,c £ {1, . . . ,n} pairwise distinct. Now choose m e {1, . . . ,n} \ {i, j, fc, I}. 
Then ip ■ {i, m} must meet each of the sets {a, b}, {a, c}, and {b, c}. Being itself a set of two elements, 
(p ■ {i, ni} must be one of the sets {a, b}, {a, c}, or {&, c}, contradicting the injectivity of (p. Therefore 
'p-{hj}'^'P-{hk}nip-{iJ}^9. 

Fix an index i e {1, . . . , n} and choose j,k £ {1, . . . , n} \ {i}. Then ip ■ {i, j} O ip ■ {i, k} is a set with 
one element, and by the above this one element must also lie in every (p ■ {i, 1} with I e {1, . . . ,n}\ {i}. 
Hence Hi^^i 'P ' {h 1} 7^ 0- This allows us to define a map a: {1, . . . , n} ^ {1, . . . , n} with 



ip-{i,j}np-{i,k} ^iD. 



(1.1) 



n 




(1.2) 



For i e {l,...,n} define Mi := {{i,j} | j G {1, . . . , ti} \ {i}}. Then (tO tells us that (p ■ Mi C M„(^i). 
Since \Mi\ — | and since ip is injective, this implies (p ■ Mi — Mg.(^iy Take i, i' E {1, . . . , n} with 

a{i) — cr{i'). Then ip ■ Mi — p ■ Mi', which implies Mi — Mi> and therefore i = i'. Thus a is injective. 
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Equation il.2\ implies that for i.j £ {1, .. .,n} distinct we can write ip ■ = {o'(i), 7i(j)} with 

7^: {1, . . . , n} \ {i} — > {1, . . . , n}. But applying (ll.2> with the roles of i and j interchanged yields 

n n 

Mj)}-n'^ =nM*),7.(j)}. 

i=l i=l 

By the injectivity of a this implies a{i) = for all i ^ j. We conclude that ■ {i, j} = {cr(i), cr(j)} 
for all i, j e {1, . . . , n} distinct. But this means that is a relabeling, as claimed □ 

Remark. For n = 4, Lemma [l~2l becomes false. An example is given hy ip E S^i^^ defined as 

^•{1,2} = {1,2}, (p. {1,3} = {1,3}, ^•{1,4} = {2,3} 
^•{2,3} = {1,4}, (^.{2, 4} = {2, 4}, ^ • {3, 4} = {3, 4}. 

This permutation satisfies ( II. Il l, but it is not a relabeling. Lemma [l~2l becomes true for 7i = 4 if we add the 
additional condition 

(^•{l,2}n^-{l,3}n(p-{l,4}7^0. (1.3) 

Do non-reconstructible point configurations exist? The answer is yes. Some examples can be found 
in Boutin and Kemper [4]. Fortunately, non-reconstructible configurations are rare. The key to this fact 
is contained in the functional relationships between the pairwise distances of a point configuration. These 
relationships are well-known from classical invariant theory. For example, a planar configuration of four 
points Pi, pj, pk, and pi satisfies 

d{,j} - - c?{j,/} -2c?{j,i} d{j^k} - d{jA} - d^k,l} 

We can also express this relationship as follows. Define the polynomial 

giU, V, W, X, Y, Z) := 2U^Z + 2UVX - 2UVY - 2UVZ - 2UXW - 2UXZ + 2UYW- 
2UYZ - 2UWZ + 2UZ^ + 2V^Y - 2VXY - 2VXW + 2VY^ - 2VYW- 

2VYZ + 2VWZ + 2X'^W - 2XYW + 2XY Z + 2XW'^ - 2XWZ. 

Then 

d{i,fc}, d{j^A;}, C?{i,i}' (^-^^ 

For simplicity, we continue to concentrate on the planar case m = 2 although other dimensions can be 
treated similarly. Recall that V denotes the set of pairs V — {{*■ j}|* 7^ j- i = 1, ■ • • i n}. The following 
theorem gives a practical test for reconstructibility of planar point configurations. 

Theorem 1.3. Letn> 5, let pi, . . . ,Pn G M.^ and let d^^jy — \\pi~pj\\^ be the square of the Euclidean dis- 
tance between Pi and Pj, for every {i, j} G V. Suppose that for each choice of indices *o>*i>*2.ii.i2.fci.fc2.^i. 
l2,mi,m2 G {1, . . . , n} such that the pairs {ig, {iq, ^2}, {^1,^2}, {^i, fe}, {^1, ^2}, {^i, '712} G V are 
distinct, we have 

g (c^{»o,»i}''^{ii,i2}''^{fci,fc2}''^{ii,i2}'C?{mi,m2},c;{»o,»2}) 7^ 0- (1-5) 
Then pi, . . . ,p„ is reconstructible from distances. 

Proof. Let qi, . . . ,qn G be a point configuration with the same distribution of distances as pi, . . . 
Write <^{j ji, = \\qi — QjW^- Then there exists a permutation ip G S/n\ of the set V such that 
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We wish to use Lemma [L2l for showing that Lp ^ is a relabehng, which will imply that ip is also a relabeling. 
Take any pairwise distinct indices i, j, fc, Z e {1, . . . , n}. Then the above equation and ( ll.4> imply 



9 {d:p-{ij},dy^.{i^k}id,^.[i^i^,d^.^j kj,d^.^j ij,d^._^k,i}) 



9 (c'lij}; '^{i,fc}i '^{ij}; C'lj.fe}; c'ij,/}; C'ifcj} 



= 0. 



It follows from the hypothesis jl.5t that (p ■ {i, j} and (p ■ {k, 1} are disjoint (otherwise they would have an 
index io in common). So for disjoint sets {i,j} and {k, 1} we have that ip ■ and ip ■ {k, 1} are also 

disjoint. This is equivalent to saying that if (p ■ {i,j} and (p ■ {fc, 1} have non-empty intersection, then the 
same is true for and {fc, I}. Take a,b,c ^ {1, . . . ,n} pairwise distinct and set {i,j} := (p~^ ■ {a, b} 
and {j,k} :— ip^^ ■ {a,c}. Then ip ■ n (p ■ {k,l} — {a, 5} n {a, c} — {a}, hence, as seen above, 

and {k, 1} have non-empty intersection. Thus the condition (II. 1> of Lemma ll.2l is satisfied for (p^^. It 
follows that (p~^, and hence also ip, is a relabeling: (p ■ {i,j} = {7r(i), Tr{j)} with tt e 5„. Now it follows 
from ProDOsition lO. ll that there exists a rigid motion (M, T) such that 

for alH S {1, . . . , n}. This completes the proof. □ 

Remark. Take indices io, ii, «2,ii,i2, ^i, ^2, ^i, ^2, ^i, ?7i2 G {1, . . . , n} as in the hypothesis of Theo- 
rem ll.3l Explicit computation shows that 

9 {d{ioM}'d{ji,j2}^d{ki,k2}^d.{i^.h}'d,{mi,m2},d,{igA^}) , 

viewed as a polynomial in variables contains the term 2d|j^ ii}'^{io-i2}- Notice that the index io occurs 

three times in this term (when writing it out as a product rather than squaring the first variable). It follows 
from Boutin and Kemper [4, Proposition 2.2(b) and Lemma 2.3] that this term does not occur in any relation- 
ship of degree 3 between the d^.j^^. In particular, g (o?{jo.ji}, o'bi ^^{fci.fca}: ^^{/i^ia}, ^{mi,™^}, ^{jo^ja}) 
is not a relationship between the It follows that there exists a dense, open subset 57 C (R-^) such 

that for all point configurations (pi, . . . e fl the hypotheses of Theorem ll.3l are met. This provides a 
new proof for the fact that "most" point configurations are reconstructible from distances, which appeared in 
greater generality in [4, Theorem 2.6]. < 

How many tests do we have to conduct for checking that the conditions in il.5i are satisfied? There are n 
choices for io, the index that is repeated. For each choice of io, there are {n — l)(n — 2) choices for ii and 
12 (since these three indices must be distinct). Having chosen io, ii, and i^, there are (2) — 2 choices for the 
set {^1,^2}, (2) ~ 4 choices for the set {fci, fc2} and so on. Altogether, we obtain 

"<"-)<"-)(0-)(0-)(0-)(0- 

— (n" - 7n^" ~ + 138n* - 83n^ - OSSn*^ + 1074n^ + 2996n* - 367271^ - 3296^^ + 3840n) 
16 

choices. 

Corollary 1.4. There exists an open and dense set Vl C (M'^)" of reconstructible n-point configurations and 
an algorithm in 0{n}^ ) steps to determine whether any (pi, . . . ,p„) G (M.^)" lies in il. 



Remark 1.5. The algorithm given by Theorem ll.3l can be generalized to M'" if n > m + 2. For each choice 
of m + 2 indices io, . . . , im+i we have the relationship 
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n 


# combinations 


CPU time in seconds 


5 


100,800 


72 


6 


2,059,200 


1,170 


7 


19,535,040 


9,920 


8 


120,556,800 


58,375 



Table 2: Time required to check for the reconstructibility of an n-point configuration. 

which can be expressed as gm (i^{io-n}' ■ ■ ■ ' ~ with an appropriate polynomial in k 

("'^^) variables. Now we obtain a generalization of Theorem ll .3l which says that if for all pairwise distinct 
choices Si,...,Sk eV with Si n Sk ^ 9 we have 

gm{dsi,---,dsj ^0, (1.6) 

then the configuration pi, . . . is reconstructible from distances. We see that there are 

n{n -l){n- 2) H ( - j) = O (n™^+3"+i) 

steps for checking the reconstructibility of pi, . . . It also follows from Boutin and Kemper [4, Propo- 
sition 2.2(b) and Lemma 2.3] that there exists a dense open subset C (M™)" where the inequalities M.61 
are all satisfied. 

2 Numerical Experiments 

A simple Matlab code (see the appendix) was used to check for the reconstructibility of some n-point con- 
figurations. In the code, we traded simplicity for speed in an attempt to make the algorithm more easily 
understandable. Even so, we were able to show that some n-point configurations were reconstructible, 
with n = 5, 6, 7 and even 8 in a reasonable time. Corresponding CPU times and number of combina- 
tions to be checked are given in Table |2] The computations were done using Matlab version 6.1 on a Sun 
(4xultraSPARC-ll, 480 MHZ). 

An important point to observe is that if a point configuration fails to satisfy one of the conditions in 
\\.5\ . it does not mean that it is not reconstructible. For example, it is not hard to show that every square is 
reconstructible (see Boutin and Kemper [4, Example 2.12]). But, as one can check, squares satisfy neither 
M.5\ nor M.2>\ . This is due to the fact that squares have repeated distances. Indeed, any planar n-point 
configuration with repeated distances will fail the reconstructibility test. (See Boutin and Kemper [4] for a 
proof of this fact and ideas on how to modify the algorithm to take care of point configurations with repeated 
distances.) Also, the point configuration given by 

PI = (0,0),p2 - (7,0),P3 - (5,-l),P4 = (3,-3),p5 - (11,2) 

does not satisfy M.5\ . even though its pairwise distances are all distinct. However, one can show that it is 
actually reconstructible. (It suffices to show that the permutations of the distances which make g equal to 
zero all violate one of the relationships that exist between the pairwise distances of five points. We checked 
this numerically.) Our test is thus not perfect. 

Observe that, when using points with small integer coordinates, the polynomial g can be evaluated exactly 
on a computer We can thus determine precisely whether such a point configuration satisfies the conditions 
of M.5\ . An interesting question is: given a planar n-point configuration with integer coordinates and lying 
inside the box [0, N] x [0, N], what are the chances that it will fail the reconstructibility test? Numerical 
experiments showed that it is quite likely, even when configurations with repeated distances are excluded. 
For iV = 3, we found that about 61% of configurations of four points whose distances are not repeated 
fail the test. (More precisely, we generated all possible pi — {xi,yi), p2 — (2^2, 2/2), Ps = (2^3, J/s), 
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P4 = (x4, 1/4) with coordinates in {0, 1,2,3} and such that either Xi < Xi+i or Xi — Xi^i and yi < j/i+i, for 
alH = 1,2,3,4. Of those 1820 four-point configurations, we found that 1636 had repeated distances while a 
total of 1748 failed the test.) For iV = 4, this percentage went down to about 30%, which is still quite high. 

It would be interesting to determine whether such high rates of failure are also observed when the coor- 
dinates of the points are not necessarily integers. But, in general, floating-point arithmetic prevents us for 
determining whether a polynomial function is exactly zero. We must thus replace the g = in conditions 
1 and 2 by |g| < e, for some e determined by the machine precision and possible noise in the measure- 
ments. However, numerical tests have shown that if the coordinates of four points are chosen randomly in 
(0, 1) (using the Matlab rand function), then the polynomial g in il.5\ rarely takes very small values. For 
example, after generating 5000 different random four-point configurations, we found that only 22 of those 
generated a g with a value less than 10"^. In another set of 5000 four-point configurations, we found only 
6 which generated a g with a value less than 10^*. In a final set of 10,000 four-point configurations, we 
found none which generated a g with a value less than 10~^. As these values are well above the maximal 
error expected with such data when evaluating g using Matlab, this implies that none of the 20,000 random 
four-point configurations we generated could possibly fail the test. 

3 The Case of Orientation Preserving Rigid Motions in the Plane 

In the previous two sections, we considered the case where the shape of an n-point configurations is defined 
hy pi, . . . ,pn £ R™ up to rigid motions. Recall that the group of rigid motions in R™, sometimes called the 
Euclidean group and denoted by E{m), is generated by rotations, translations and reflections in R™. How- 
ever, in certain circumstances, it may be desirable to be able to determine whether two point configurations 
are equivalent up to strictly orientation preserving rigid motions. The group of orientation preserving rigid 
motions, sometimes called the special Euclidean group and denoted by SE{m), is the one that is generated 
by rotations and translations in R™. 

For simplicity, we again restrict ourselves to the planar case m = 2. Given a planar point configuration 
pi, . . . ,Pn S R^, we would like to be able to determine whether any other planar n-point configuration 
qi, . . . ,qn is the same as pi , . . . , p„ up to a rotation and a translation? Given any qi , qj , q^ in the plane, 
denote by aq^^q- the signed area of the parallelogram spanned by qi — qj and q^ — qj, so 

= det(gj -qk,qj -qk)- 

Since signed areas are unchanged under rotations and translations, the function / : 
defined by 

("^91.92,93 ~ "^gi, 92, 94) ("^91 '92, 93 ^ '3^9l'92,94 +2agj^g3^g^) 

('^91,92,93 ~ 2agj^g2,g4 + , gg , ^4 ) (2ag j ^ ^ gg — 0,qi,q2,qi + ^91,93,94) (3-1) 

is invariant under the action of SE{2). Moreover, one can check that it is also invariant under a relabeling 
of the four points qi,q2, q3, qi- However, it is not invariant under rigid motions in general. Indeed, any 
transformation which is a rigid motion but does not preserve the orientation will transform / into — /. 

Given an n-point configuration qi, . . . ,qn with n > 4, we can evaluate / on all possible subsets of four 
points of {qi, . . . , qn}- We consider the distribution of the value of these /'s, i.e. the distribution of the 

liiMMM = I{qii,qi2^qi3:qii)^ for aii ii < 12 < ^3 < «4 e {1, • • ■ 

Proposition 3.1. Let n > 4 and Zef pi, . . . ,p„ g be an n-point configuration which is reconstructible 
from distances. Assume that the distribution of the I's of this point configuration is not a symmetric function 
(i.e. that the distribution of the I's is not the same as the distribution of the —I's.) Let qi, G be 

another n-point configuration. Then both the distribution of the distances and the distribution of the I 's of 
the two point configurations are the same if and only if there exists a rotation and a translation which maps 
one point configuration onto the other 
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Proof. Observe that, in addition to being invariant under rotations and translations of the points, the dis- 
tribution of the value of the /'s is also independent of the labeling of the points. The same holds for the 
distribution of pairwise distances. So if two n-point configurations are the same up to a rotation, a transla- 
tion and a relabeling, then the distribution of the /'s and the distribution of the distances are the same for 
both. Thus the //is clear. 

Now assume that the distribution of the distances and the distribution of the J's are the same for both point 
configuration. Since pi, . . . ,pn is, by hypothesis, reconstructible, this implies that there exists a rigid motion 
(M, T) and a relabeling tt : {1, . . . ,n} ^ {!,...,«} such that Mpi + T — qT^(i), for alH — 1, . . . , n. 
If (M,T) is not in SE{2), then it maps each I{Pn,Pi2,Pis,Pii) to -I{pi^,pi^,Pi^,piJ). But this is a 
contradiction, since the distribution of the I's is not symmetric. Thus g is in SE{2). This shows the only 
if. □ 

Remark 3.2. One can actually show that if pi, . . . ,p4 is equivalent to gi , . . . , 54 up to a rigid motion, 
then pi, . . . ,p4 is equivalent to gi, . . . , 94 up to a rotation and a translation if and only if I{pi, . . . ,^4) = 
I{qi, . . . , (74). (Indeed, / is one of the two fundamental invariants of the action of SE{2) x 5*4 on ^ 



p2 



X which we obtained using the invariant theory package in Magma [3]. By construction, these two 
invariants thus distinguish the orbits of SE{2) x ^3. The other invariant is actually unchanged under the 
action of the full Euclidean group E{2) and so / alone distinguishes the orbits of SE{2) within the orbits of 

E{2). 



4 The Case of Rotations, Translations and Scalings 

In certain circumstances, it may also be desirable to be able to determine whether two point configurations 
are the same up to a rigid motion and a scaling. This can be done using a simple variation of the previous 
approach. Given a distribution of distances = \\pi — PjlP}, let dmax be the largest distance 

draax = max{d{i | {i , j} e V}, 

which can be assumed to be non-zero since otherwise all points coincide. We can consider the distribution 
of the rescaled distances {-;^^-^}{i.j]e'p- I" addition to being invariant under rigid motions and relabeling, 
the distribution of the rescaled distances is also invariant under a scaling of the points 

Pi ^ Xpi, for every i = 1, . . . , n. 

for any A G M^o- 

Proposition 4.1. Let n > m + 2. There exists an open, dense subset i7 o/(M™)" such that if an n-point con- 
figuration pi, . . . ,pn is such that (pi, . . . ,pn) G fl, thenpi, . . . ,pn is uniquely determined, up to rotations, 
translations, refiections, scalings and relabeling of the points, by the distribution of its rescaled pairwise 

distances { rf''"'^ }{»j}g?- Moreover, there is an algorithm in 0{n t""*^ ) steps to determine whether 

(pi, ■ ■ ■ ,Pn) e n. 

Proof Let pi, . . . , pn G K"' be an ?i-point configuration which is reconstructible from distances and whose 
pairwise distances are not all zero. Observe that if gi, . . . , g„ S is another n-point configuration, then 
the distributions of the rescaled distances of both point configurations are the same if and only if there exists 
a rigid motion followed by a scaling which maps one point configuration onto the other. The claim is thus a 
direct corollary of Theorem 2.6 from Boutin and Kemper [4] and of Remark [T31 □ 
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Appendix 

% This is a simple Matlab function that determines whether 
% the pairwise distances between the points of the planar n-point configuration 
% defined by the n columns of a 2-by-n matrix p satisfy the conditions of (II. 5> 
% with replaced by some (small) number e. 

function g=evaluate_g(p,e) 
% p is a 2-by-n matrix. 

% e should be chosen depending on the machine precision. 
% This function returns if 5 < e for all sets of pairs 
% of the conditions in ll.Sl and 1 otherwise. 
n=length(p); 

% Compute the squares of the distances d. 
fori=l:n-l, for j=i+l:n 

d(ij)=sum((p(i,:)-pa,:)).^2); 
end, end 
d(n,n)=0; 
d=d+d'; 

for iO=l:n, for il=l:n, if ir=iO 
xl2=d(i0,il); 

for i2=l:n, if i2~=il & i2~=i0 
x34=d(i0,i2); 

for jl=l:n-l, for j2=jl+l:n 

if (jr=iO I j2~=il) & (jr=il I j2~=i0) & (jr=iO | j2~=i2) & (jr=i2 | j2~=i0) 
xl3=d(jl,j2); 

forkl=l:n-l,fork2=kl+l:n; 

if (kr=iO I k2~=il) & (kr=il I k2~=i0) & (kr=iO I k2~=i2) & (kr=i2 | k2~=i0) & . . . 

(kr=ji I k2~=j2) & (kr=j2 1 k2~=ji) 

xl4=d(kl,k2); 
forll=l:n-l,for 12=ll+l:n 

if (ir=io I i2~=ii) & (ir=ii I i2~=i0) & (ir=io | i2~=i2) & (ir=i2 1 i2~=i0) & . . . 

(ir=jl I 12~=j2) & (ir=j2 I 12~=jl) & (ir=kl 1 12~=k2) & (ir=k2 I 12~=kl) 
x23=d(ll,12); 

for ml=l:n-l, for m2=ml+l:n 

if (mr=iO I m2~=il) & (mr=il | m2~=i0) & (mr=iO | m2~=i2) & . . . 

(mr=i2 I m2~=i0) & (mr=jl | m2~=j2) & (mr=j2 | m2~=jl) &. . . 

(mr=kl I m2~=k2) & (mr=k2 | m2~=kl) & (mr=ll | m2~=12) & (mr=12 | m2~=ll) 

x24=d(ml,m2); 

mll=-2*xl4; 

ml2=xl2-xl4-x24; 

ml3=xl3-xl4-x34; 

m22=-2*x24; 

m23=x23-x24-x34; 

m33=-2*x34; 

Mu=mll*m22*m33-mll*m23"2-ml2"2*m33-ml2*m23*ml3+ml3*ml2*m23-m22*ml3"2; 
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ifMu==0 
g=0; 



return 



end, end, end, end, end, end, end, end, end, end, end, end, 
end, end, end, end, end, end 
g=i; 
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